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THE TURBULENT FLOW IN DIFFUSERS OF SMALL DIVERGENCE ANGLE 1 
By G. A. Gourzhienko 

SUMMARY 


The turbulent flow in a conical diffuser represents the type of 
turbulent boundary layer with positive longitudinal pressure gradient. 
In contrast to the boundary layer problem, however, it is not neces- 
sary that the pressure distribution along the limits of the boundary 
layer (along the axis of the diffuser) be given, since this distribu- 
tion can be obtained from the computation. This circumstance, together 
with the greater simplicity of the problem as a whole, provides a use- 
ful basis for the study of the extension of the results of semiempiri— 
cal theories to the case of motion with a positive pressure gradient. 

In the first part of the paper, formulas are derived for the com- 
putation of the velocity and pressure distributions in the turbulent 
flow along, and at right angles to, the axis of a diffuser of small 
cone angle. The problem is solved on the basis of the following 
assumptions: 

1. The motion is assumed to take place along straight lines inter- 
secting at the vertex of the diffuser cone . 

2. The normal components of the turbulent stress tensor are as- 
sumed isotropic. Their gradients along the diffuser are neglected by 
comparison with the gradient of the static pressure. 

3. In the equations of motion the hypothesis of the mixing length 
(Prandtl formula) is applied, and it is assumed that the curve of de- 
pendence ef the nondimensional mixing length on the distance from the 
wall is absolute. 

4. In determining the shape of this curve none of the existing 
turbulence theories is taken as a basis, the assumption, common to all 
turbulence theories, only being made that the first derivative of the 
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mixing length with respect to the distance from the wall is the nondi- 
mensi onal universal turbulence constant x . The shape of the mixing 
length curve is chosen as a cubic parabola from the boundary condi tions 
applicable to it. 

5. To- obtain the solution of the equation of motion in finite form 
there is applied the step-by— step interpolation of the values of the 
inertia integral, the applicability of which is proved. 

6. At the walls of the diffuser the existence of a laminar sub- 
layer is assumed, the thickness of which follows the known Karman law. 
It is shown that the assumption of radial motion in this sublayer is 
fundamentally inapplicable. In obtaining the velocity distribution in 
the sublayer, the assumption is made of the continuity of the curve of 
friction distribution in passing through the boundary of the sublayer. 

7. The resistance formula obtained for the diffuser is found to be 
absolutely identical with that for the pipe . 

In the second part of the paper tests are described on the meas- 
urement of the velocity and pressure distribution in two conical dif- 
fusers with angles of 1° and 2°, and a detailed comparison is made 
between the experimental and the theoretical results. It is found that: 

1. The assumption of the radial character of the flow is satisfied 
only with a certain, though large, degree of approximation as should be 
the case if viscous friction exists simultaneously with turbulent fric- 
tion. 


2. The obtained formula for the velocity distribution agrees well 
with the experimental results. 

3. The increase computed (on the basis of the derived resistance 
formula with the values of the universal constants taken from tests on 
pipes) in static pressure along the diffuser deviates little from the 
test results. 

In the third part of the paper a first, very approximate, at- 
tempt is made at estimating semi empiric ally the deviation of the true 
motion from the radial pattern assumed. The analysis is based on the 
assumption that the true motion may to a first approximation be as- 
sumed radial but emanating from another fictitious source . The 
latter, for the region near the diffuser axis, is computed by two 
different methods. The good agreement of the results of these com- 
putations shows that the assumption is justified. The assumption of 
the displaced source in computing the characteristics of the growth 
of the static pressure along the diffuser gives complete agreement 
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with the experiment for the 1° diffuser, while for the 2° diffuser 
the disagreement is greater than without the assumption. The result 
for the 2° diffuser evidently is explained hy the fact that the gradi- 
ent'' of the normal component of the turbulent stresses was not taken 
into account together with the pressure gradient. This correction, 
which is negligible for the 1 diffuser, may have an appreciable effect 
in the case of- the 2° diffuser. 


INTRODUCTION 


The chief experimental source on which at the present time the so- 
called semi empirical theories of turbulence are based is the fully 
developed turbulent flow between two parallel planes or in a straight 
cyli ndm cal pipe . The reason for such an exclusive role played by 
these two types of flow lies in the circumstance that these are the 
simplest types of flow as regards their kinetic and dynamic relations. 
In the first place, both for the case of flow between parallel walls 
and for the circular pipe there is no need to consider the change in 
the mean velocity and friction profiles along the axis of the flow 
since neither the velocities nor the frictional stresses along the flow 
direction change •' n value. In the second place, because of the absence 
of inertia forces the profile of the shear stresses transverse to the 
flow direction is found bo be linear’. This fact is not a consequence 
of an, hypothesis regarding the turbulence but follows from the funda- 
mental equation of notion (Reynolds), finally, for these sample cases 
the equation of motion permits the experimental computation of the 
shear stress at the wall by measuring the drop in static pressure along 
the flow, as may be done vith very great accuracy. The latter circum- 
stance very greatly simplifies the experimental confirmation of the 
theoretical resistance laws. 

The careful experimental investigation of the above two cases of 
turbulent flow has led to a completely satisfactory application for 
practical purposes of the semiempirical turbulence theories (Prandtl, 
Karman, Mattioli) . The empirical nature of these theories lies, as is 
known, in the presence of the experimental constants X (the absolute 
turbulence constant) and a (the nondime ns lonal thickness of bhe lam- 
inar sublayer at the wall) obtained in evaluating the test data on the 
relation between the friction at the wall and the Reynolds number from 
the formulas obtained on the basis of these theories. 

At the present time it may be confidently asserted that f»r the 
above-mentioned simplest cases of flow the constancy and. absoluteness 
of these constants are facts that have been repeatedly verified. How- 
ever, it still remains very uncertain to what extent the absoluteness 
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of these constants is maintained in passing to other cases of motion, 
in particular to motion, not with a negative static pressure gradient 
(as in the pipe and channel) , hut with a positive gradient — that is, 
to those cases characterized h y the existence of inertia forces. As 
is known, in these cases both the velocity and the frictional stress 
profiles undergo radical changes. Such a case is that, for example, 
of the turbulent motion in the boundary layer at the upper surface of 
a wing. Notwithstanding the undoubted urgency of the solution of the 
problem of the wing boundary layer this type of motion appears very un- 
suitable for the purposes of generalizing the existing semi empirical 
turbulence theories. Thus, the frictional stress along the wing con- 
tour varies according to a law which is connected with the velocity 
distribution law by a very complicated integro-differential condition 
(Karman) . 

In addition to the fact that this condition on evaluating the test 
data requires the carrying out of graphical differentiation of the ex- 
perimental curves, a procedure which introduces an element of arbitrar- 
iness in the case of the flow about a wing, there, strictly speaking, 
does not exist a determinate problem sinc.e the static pressure distri- 
bution over the wing profile must in all cases be obtained from exper- 
iment. Thus the "turbulent" character of the phenomenon in the 
boundary layer is very mu,ch complicated by accessory circumstances of 
the external problem. For this reason it was considered desirable to 
obtain a type of flow which, while possessing all the properties asso- 
ciated with a variable positive pressure gradient, was most free from 
various external complicating circumstances. 

Such a type of motion, that is, a somewhat exaggerated model of 
the boundary layer of a wing, may be represented by the steady flow in 
a straight— walled or conical diffuser. Since m this case the entire 
region within the diffuser is filled with the "boundary layer," the 
chief difficulty of the external problem - namely, the Incompleteness 
of the equations of motion as regards the external conditions — drops 
out. The equations obtained are determinate both for the velocity and the 
pressure distributions. Moreover, by introducing a certain assumption, 
for small divergence angles of the diffuser there is the possibility of 
greatly simplifying the investigation of the character of the change in 
the values of the velocities, friction, and pressures along the flow by 
entirely avoiding the operations associated with graphical differentia- 
tion. By 'investigating the possibility of generalizing the semiempirical 
theory to the case of the diffuser it is possible then to proceed on a 
surer basis to the boundary layer study. In selecting the conical dif- 
fuser for the present investigation the following facts were considered 
because of convenience in testing: 
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1. To obtain steady turbulence in the diffuser, it is necessary to 
have an initial inlet length ahead of it. It was convenient for this 
purpose to use the cylindrical pipe, already investigated in a previous 
paper (reference 1), which may be considered as a "diffuser of zero 
divergence . " 

2. In investigating a straight-called diffuser with an initial in- 
let length, it is necessary to make the distance between its lateral 
parallel walls as large as possible in order to avoid the effect of in- 
creasing boundary layer on the flow in the axial plane . A large dis- 
tance between the lateral walls would give such a large diffuser cross 
section, however, that the air intake apparatus at disposal would 

be unable to produce a flow with sufficiently large Reynolds number. 

The present paper is divided into three parts. In the first part 
an attempt is made to give a theoretical analysis of the turbulent flow 
in a conical diffuser with small divergence angle. In the second part 
the procedure is described, the results of the tests conducted on coni- 
cal diffusers with cone angles of 1° and 2° are presented, and an 
exhaustive comparison of the experimental results with the theory pre- 
sented in the first part is- given. In the third part an approximate 
method is given for estimating the deviation of the true motion from 
the radial pattern assumed in the first part. 


I . THEORETICAL ANALYSIS OF THE TURBULENT FLOW IN A CONICAL DIFFUSER 
1. Fundamental Equations of Motion 


It will be assumed that a fully turbulent flow enters from a 
straight cylindrical pipe of radius R into a conical diffuser with 
angle between the axis and the generator (fig, 1). To investigate 

this case, it is convenient to assume a spherical system of coordinates 
with the pole at the vertex of the con© 0 and with the polar axis 
directed along the axis of the cone. Thus the coordinates of any point 
M within the diffuser will be the distance from the pole r and two 
angles, namely, -3 between the straight. line connecting the point with 
the pole and the polar axis, and cp between the perpendicular dropped 
from the point M on the polar axis and any fixed plane containing the 
polar axis. 

The hydrodynamic equations (Euler) in the chosen coordinates have, 
as is known, for the steady motion the following form; 
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l > ( 1 ) 


1 Sj> 
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and the equation of continuity 


Sv r 1 

Sv.,3 ^ i ' 


Sr r 

S8 r sin 8 
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^8 


+ — cot 8=0' 
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( 2 ) 


where v^ is the component of the velocity directed along the straight 

line passing through the pole and v.» and v-p are the velocity com- 
ponents directed along the tangents to the arcs measuring the corre- 
sponding angles; p is the static pressure in the flew, and p the 
density of the fluid. From t'-'p pro. '“ding equations is set up the equa- 
tion of Reynolds for which purpose, as usual, the velocity and pressure 
are broken down into a mean val^e r th respect to time and a fluctua- 
tion about the mean: 


v S3 V + V 
v p v I* I* 




v cp “ v q> + v cp 


(:•) 


v 8 = v 8 + v 8 


P - P + P 1 


Equations (1) with the aid of the continuity equation may be trans- 
formed into the following; 
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r sin a ocp r 


2 


cot a 


_ 1 dp 

pr sin a cep 

By substituting (3), averaging with respect to time and remember- 
ing that according to the averaging laws: 


v r v a = ya + v r’ v a ' 

v* = 0, etc. 

the three equations of motion are obtained: 
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and the continuity equation 
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T rcp 

T <pd 


P ^ 


P v r v tp 


p v* 


> ' . 


are the tangential and 



t 3-3 


.2 

p 




i 


1 2 

T = — p V 

W 9 


the normal components of the stress tensor. 

Since the Reynolds eo nations were derived from the Euler equations, 
that is, equations that do not „.ake the viscosity into account, it is 
safe to assume that the tangent al components of the stress tensor 

T r$ ' T rcp> and V 

represent the sum of the viscous and turbulence parts of the stresses, 
because the viscous stresses, by analogy with the turbulence stresses, 
are also a consequence of the averaging with respect to time of the ac- 
tual molecular motion. 

If the analogous components of the viscous stresses are substituted 
(expressed in spherical coordinates in terms of the velocity gradients) 
in the obtained equations in place of the tangential components of the 
turbulence stresses, and if isotropy of the normal viscous stresses is 
assumed, the wo 11— known equations of Navier-Stokes are obtained. 

As is known, the system of Reynolds equations is not determinate 
since, for determining the ten unknown functions. 
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r q> } r v-3’ T rcp , 


T rr' 




cpcp 


only four equations are available. Some simplifying assumptions now 
will be introduced. It will be safe to assume that the average motion 
in the diffuser is one with axial symmetry; that is, the velocity com- 
ponent v^ and all derivatives with respect to <p are equal to zero. 

This assumption immediately removes the third equation from considera- 
tion and simplifies the others. Then, introduce the less obvious 
assumption that the mean motion is along the straight lines passing 
through the pole. This, of course, is a very strong assumption, approx- 
imately satisfied evidently only for small divergence angles of the 
diffuser, and, no doubt, requiring experimental confirmation (which will 
be presented later) . By this assumption it may be considered that 
v^ = 0. In this case it follows that: 


2v 
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p Sr 


1 | 1 Sx rr 
p r or 


1 - T rfl _ r jn _ + + Ifi d 
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Introduce the assumption of isotropy of the normal components of 
the turbulent stresses, that is, that the equations are approximately 
satisfied 


T 

rr 
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Then there is obtained from the first equation, letting for simplicity 
v r = u and T r3 = T 


_ Su 2u" 

2u — + 
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where p — T rr represents some total pressure which at the diffuser 
wall where the fluctuations vanish, is equal to the static. Setting 
P “ T rr = Pi> and assuming in view of the smallness of the angles 
sin d d gives: 

2u^+^=~i ~ + JL Mlil (k) 

Sr r p Sr S-a 

The second equation containing the gradient of p transverse to 
the difftvser is excluded entirely from consideration. This equation is 
obtained by taking the sum of the projections of all forces acting on 
the fluid in the direction perpendicular to the mean radial motion. 
Therefore in introducing the 'assumption of the radial character of the 
motion, this equatioii suffers to a considerably greater degree than the 
first, obtained from considering the projections of all the forces in 
the direction of the initiai motion. It must be supposed that the 
neglected terms in the second equation, due to the assumption of radial 
motion, which contain the component v^ are of the same order of mag- 
nitude as the remaining terms containing — and the difference 

dd 

T cp 9 — T ■Q’Q • Thus the retention of the second equation must be considered 
as unsuitable and even harmful. In what follows it will be seen tha,t 
the remaining equations are entirely sufficient for solving the problem 
for the assumptions made. 

The continuity equation under the assumption of axial symmetry and 
radial character of the mean motion becomes: 



dr r 


This equation may be easily Integrated with respect to r and its gen- 
eral integral will be: 


= Mil 

r 2 


(5) 


where t 


is an arbitrary function of d only. Expression (5) shows 
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that under these assumptions- similar velocity profiles should he obtained 
over the entire length of the diffuser and that along each radial line of 
flow the velocity varies in the same manner as for the motion of an ideal 
fluid - that is, inversely proportional to the square of the distance 
from the source. 

Substituting equation (5) in the initial equation (4) results in 
the following equation after combining similar ' terms and multiplying by 
r 5 : 

2t s a ~ r s _ JL d(-rfl) 

p dr , pd d-fl 


From the equation obtained two very important conclusions can be drawn 
on the change in pressure and intensity of friction along the diffuser. 
Since the left side of the equation does not depend on r, it must be 
assumed that the right side likewise is independent of r. This can be 

-the case only if cpi/dr is inversely proportional to r 5 and r is 

4 

inversely proportional to r . 

By setting 


p 1 &Pl 5 

G = — r 


where G and f according to what was said above do not depend on r, 
the above equation of motion may be written as the usual one: 



cL(fO) 

diS 


(7) 


The first expression in (6) may easily be integrated with respect to r. 
There is obtained 
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where. p lQ is.th^ absolute value _ of the generalized- pressure at any 
poipt of the flow having’ the coordinate ' r^. The’ expression. (8) giving 

the law of variation of this pressure along the diffuser, together with 
the law of drop in velocity (5 ), is very suitable for the experimental 
checking of the fundamental assumption of radial flow. Use will be made 
of this expression later. 

It was assumed that .the,, pressure. p x is the difference between the 
static pressure and the normal turbulent stress:- 


Pi 9 P T j. r 


Thus the magnitude G-, strictly speaking, is determined as 

G - S 

p V Sr 


dr 


rr 


dr 




It is assumed that the gradient along r of t is small by compari- 
son with the gradient of the pressure p. This makes it possible to 
consider 



In studying the boundary layer the assumption is usually made that 
the static pressure does not vary over the thickness of the bpundary 
layer. This is excellently confirmed by experiment. Applying an analogous 
assumption to thiB case, let 



This at once leads to the result that G = constant. Multiplying equa- 
tion (7) by d and integrating from -8 = 0 (axis of diffuser) to a 
variable -9 gives 



t 2 ddd = 



2 
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Transform to ncndimensional variables fol* which purpose the angle is re- 
ferred. to its maximum value -d or lotting - 


= l 


The velocity u will be referred to the "dynamic velocity of friction" 


at the wall v = / t /p where t is the frictional intensity at 


the wall , and then u/v _ = rp. Since u = t/r 2 and 


V = /T /p » / f /r^p = / x 

y o' •'o' 2 / o 


/p 


t = bep is obtained where b = y'fjp. Substituting in the last ex- 
pression yields 


2 / cp 2 |dfi = ~ i- H; 

' o 2b 2 «o pt)‘ 


' , (9) 


On substituting in the second term of the right side * - / V p > 

there as obtained the generalized law of the distribution of the fric- 
tional stress transverse to the diffuser: 


f_ = t__ _ 1 
f 


r " 7 I 

r\ — v 


P f - 


cf^fdq - 


*o<H* 

2b 2 


( 10 ) 


Consider this expression. For s = 0 on the diffuser axis there 
obtains, on applying the rule of L’Hospitai, 

, t 


2d •/ 
o J 
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© 2 S;dS 


= lim 
E o 


qj^E X 2d n 

= lim — = 0 

t o 1 
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At the diffuser wall for 6=1, it is seen thdt f = f and 

’ o 


20 o f q> a 6de - 


j£ 

2b 2 


= 1 


(ID 


which gives the relation between the integrals of the inertia forces, 
the frictional forces, and the forces arising from the drop in static 
pressure . 


In order to estimate the shape of the curve of friction distribu- 
tion across the diffuser for various signs of the pressure gradient, the 

a / F ' \ 

derivative — - ( — ) will be found for ? = 1, that is, at the dif— 
5? Vf y 

o 

fuser wall. Differentiating (10) with respect to I, gives: \ 


d_ 

dt 


= d r 


2cp e 


G 

2b 2 


E 

n 


-bJ 


<p 2 ldl 


( 12 ) 


Substituting 5 = 1 and bearing in mind that at the wail cp = 0 yields 


d 

L d?- 


\ f D J - |=i ~ \ b 2 + / 


For not too small divergence angles of the diffuser, the sign of the 
magnitude in parenthesis usually is entirely determined by the first 
term. Thus, the sign of the derivative of the curve of friction dis- 
tribution will be negative for G > 0; that is, for positive pressure 

gradient and the curve of friction distribution f,/f 0 = f( p ) will have 
the form shown on figure 2. As approaches zero, that is, in 

passing from the diffuser to the straight pipe, there is obtained from 
(10) bearing in mind (11): 


T 


T 


o 
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that is, a linear distribution of the frictional stress. This straight 
line, corresponding to the negative pressure gradient, is also shown in 
figure 2. 

It is natural to assume that the method explained is suitable also 
for considering the motion in a converging pipe where i9 Q < 0. In this 

case the point of intersection of the radial lines of flow lies on the 
other side of the section considered. Since r in this case will oe 
reckoned opposite the flow direction, it follows that 


dp 

Sr 


> 0 


although the pressure drops 
pipe on the basis of (12), 

(fig- 2). 


along the flow. Thus for the converging 

d / f N - 

— ( : — j will ce positive for 5 = 1 
d ‘ Vf o' 


It is not difficult to conjecture, on the basis of the above dis- 
cussion, that a diffuser may be imagined for which there is no change 
in pressure along the flow; that is, G = 0. The divergence angle of 
such a diffuser will be determined on the basis of (li) ae x 




o 


(p=ccnst) 



(13) 


K. E. Fedyaevsky, in his paper on the boundary layer of a wing 
(reference 2), has shewn that for the externai problem in the absence 
of a pressure gradient along the surface of the body (flat plate) the 
derivative of the curve of friction distribution in the direction per- 
pendicular to the surface at the wall is equal to zero. It is inter- 
esting that for the diffuser this assumption is not f ovoid to be correct. 
According to expression (12) the derivative of the friction distribution 
at the wall in the absence of a pressure gradient along the diffuser — 
that is, for G = 0 - is equal to 


1 Such motion represents the boundary layer of an infinitely thin 
flat plate set at zero angle of attack. 
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f =1 

p=const 


= - 1 


The corresponding curve is shown on figure 2. 

Consider what the shape of the curve of friction distribution will 
he near the diffuser axis, that is, for 1=0. The first derivative 
of the curve of friction distribution for £ ~ 0 is obtained from (12) 
on substituting 1=0. Evaluating the indeterminate expression 


I 

lim / cp 2 £ d? 

?-»o J 


gives: 


I 

lim j cp 2 £d* = 2 11 m = q> £ - 

l-o 6 2 { i-o 2 ? 


On substituting - : n (11) there is obtained 






m 


2 b 2 J 


where <p m is the value of the nondimensional velocity u/v # on the 

diffuser axis. Thus, the value of the derivative of the curve of fric-r 
tion distribution on the diffuser axis and, of course, also in the 
general case is not equal to zero. In the boundary layer, according 
to the investigation of Eeayaevsky, this derivative Is equal to zero. 
This sums up the conclusions which may be drawn from the incomplete" 
equations of motion on the basis of the assumed hypotheses. 
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2. Completion of Equation of Motion 

In order to obtain the profile of the velocity distribution across 
the diffuser and the resistance law, it is necessary to render the fun- 
damental equation (9) determinate by connecting in some way the general- 
ized frictional stress f = Tr 4 with the remaining variables. This, as 
usual, is attained by introducing some semiempirical assumption with re- 
gard to the turbulence. 

It is possible, of course, to take any one of the semiempirical 
turbulence theories available at the present time (Prandtl, Karman, 
Mattioli) and generalize it to the case of the friction in a diffuser. 
This generalization to a first approximat j on (evidently quite satis— 
tory for small divergence angles) may be brought about by simply passing 
from cylindrical coordinates to polar. At first this method was followed 
in an attempt to generalize the theory of Mattioli. The unusual, com- 
plexity, however, in the solution of the differential equation obtained, 
the impossib:' 11 ty of obtaining the velocity profile in a finite form 
suitable for practical application, 1 and so forth, made it necessary to 
drop this method and seek another. 

Still another consideration led to this resolve. All empirical 
theories without exception are^built on a very shaky physical basis. 
Neither the assumption of Karman on the similarity of the fields of 
velocity fluctuations nor the assumption of Mattioli on the transport 
of momentum can at the present time be supported by any but those 
authors. It is very significant that at the present time papers devoted 
to the analysis and improvement of these theories no longer appear. All 
this indicates that interest in these theories has dropped sharply and 
that they are now in the passing stage. On the basis of what has been 
said it may be asked whether it is worth the effort, by overcoming the 
huge computational difficulties, to generalize any of the theories men- 
tioned to the case of flow in a diffuser. Even after having overcome 
all the computational difficulties and having obtained excellent agree- 
ment of the results of computation with experiment, no progress in 
learning the mechanism of the turbulent motion will have been made, for 
it would be very difficult to establish on what grounds the good re- 
sults were obtained, whether on account of the correctness of the 
generalized semiempirical theory or on account of the fact that the 
particular nature of the flow in the diffuser was already sufficiently 
taken into account by the incomplete equation of motion. It would, 
however, be a mistake to think that the semiempirical theories of 


x As is known, even for the case of flow in a pipe, the theory of 
Mattioli leads to incomputable analytical quadratures. 
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turbulence which have been mentioned have contributed absolutely nothing 
to the understanding of the physical picture of the phenomenon. It is 
sufficient to recall that all these theories lead to the same resistance 
law for the straight pipe: 


1 


/ 



x y 2 


(no 


(where 




B is the radius of the pipe, X and 


A are universal constants), which is excellently confirmed by numerous 
tests. This leads to the supposition that all these theories contain 
some common element that correctly represents the phenomenon. In order 
to explain what constitutes this common element, it is recalled that 
Prandtl connected the frictional shear with the derivative of the mean 
velocity by means of the relation 


2 

t = — p u’v' = — pZ 2 ^ (15 


where Z, a linear magnitude denoted by Prandtl as the "mixing length" 
(Mischungsweg) , is a function of the distance y from the wall of the 
pipe or channel. From its meanipg Z should become zero at the wall 
where no mixing can take place . Prandtl made the very simple and ele- 
gant assumption that follows from the possibility of developing the 
function Z = Z(y) into a series at the wall; namely, he assumed that 
to a first approximation 

Z = Xy (16) 


where X, a universal turbulence constant, is the first derivative of 
the mixing length with respect to v at the wall; that is, for y = 0 
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As is known, the assumption (15) led Prandtl to the "universal velocity 
profile" aM' to the resistance law (lk) . With the aid of the latter 
the magnitude of X was found and its absolute character established. 

Since from other theories (Karman, Mattioli) the same formula was 
obtained, it is to be expected that for these theories condition (17) 
is likewise satisfied - which statement will be proved. 

✓ / 

As is known, Karman obtains a relation connecting the mixing length 
with the derivatives of the velocity in the form: 


l 


= X 


u* 


u 


■ XE 


£l 

cp" 


and the velocity distribution in the form: 

<P m ~ 9'= - --[log(l - / t) + / t] 

X * 

By finding the derivatives cp 1 and <p" and substituting in the expres- 
sion for Z there is obtained 


l = 2XB(/ 5 - ?) 


Differentiating with respect to t gives 


dZ_, 

dl 




) 


By substituting t = 1 and remembering that df = d(l - y/E) = - dy/E 
there is obtained 



= X 


as was to be proved. 
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To Verify expression (17) according to the theory of Mattioli , 
equation ( 15 ) is solved for Z and differentiated with respect to y , 
where 



Introducing the variables cp and | fields: 

' di _ _ f 1 _ d_ ( 1 ^ yj "! 

dy - cp* S\ di V <pi J J J 


Ab is known, the theory of Mattioli for large Reynolds numbers leads to 
the following relations: 


(a) 


cp* 


X© 

ce , 


where 


© - cp — p P = constant 

£ c = constant 


(h) 


1 di _ _ 1 * 

I d© cp* + 6 


Making use of these relations gives, from the expression for dl/dy. 


dZ 

dy 


l 

/fi 


X 


'v 1 


© 


\ 

\ 


ce 


X9 j 


_1 

X© 

ce 


Since for the case where the viscosity is neglected (large Reynolds 
numbers) , the velocity at the wall approaches , 
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l'im © — » 

£-»o 


^ im ~T® = 0 

0— "5° ce 



a; X 


as was to be proved. 

It is thus evident that the fundamental element common to the fore- 
going theories is the fact that the first derivative of the curve of 
the mixing length against the distance from the wall has at the wall, 
or more accurately, at the edge of the laminar layer, a constant abso- 
lute value X independent of any variables. 

Since the resistance formula, which is a consequence of this 
assumption, is in excellent agreement with the experiment, it may be 
considered that also for the diffuser and, in general, for any turbulent 
flow near the wall 


.( 


dZ_ 

ay 


\ 

) 

y-o 


= X = constant 


Moreover, numerous experimental investigators (Conch, Nikuradse, 
Frietsche) have established two further significant fac’ts: (1) the 


curve 


B 


^ y 




f ( - 
B 


is found to be almost absolute near the wall for 


the most varied cases of flow and (2) on the axis of the channel or pipe 


or on the edge of the boundary layer 


E 


constant «« 0.l4. 


On the basis of all that has been said, it is not necessary in the 
present paper to make use of any developed hypothesis of turbulence ex- 
cept the formula of Prandtl (lb) but a method will be used that permits 
obtaining the simplest and clearest results. 
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Up to now it has teen understood that t represents", the total * 
frictional stress: * 


T = T lam + T turb 


As earlier investigations (references 1, 3, and 4) have shown, it may be 
considered that for practical Reynolds numbers the effect of the viscous 
friction is negligible in comparison with the turbulent friction, for 
this reason, in what follows, by r will be understood only its turbu- 
lent part. 

The formula of Prandtl will be reduced to a form independent of r. 
On assuming, according to the foregoing discussion, the ratio if R 
(where R is the local radius of the diffuser cross section, R = r-9 n ) 

independent of r, to be a function of -9 , from equation (15 ) there is 
obtained 



It is seen that on the basis of the assumption of the absoluteness of the 

curve — = f(5) it was found from the formula of Prandtl that the tur— 

R 

bulent part of the frictional stress, in the same manner £|s the total 
fractional stress, should be inversely proportional to r . 

The formula of Praiidtl is reduced to' nondimens I cnal form by letting 
t s= bcp; 



. d<p 
where cp* s — 
d! 


Substituting the expression obtained in the fundamen- 


tal equation of motion in the form of equation (9)- gives 



(18) 
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It ▼ill fee nsceaeary to obtain 
function 



in the form of a power 


l .n 
— = as + m 
E 


(19) 


where the unknown constants a, m, and n are obtained from the fore- 
going consi derations: 

(a) for 1=1 (at the wall) -• = 0 

E 

(b) for £ = 1 (at the wall) — * — — ) = X 

dy as \ E / 

(c) for 1-0 (at the axis) — = 0.14 

E 

By neikiijg use of the first condition there is obtained 

a + m = 0 (I) 

Differentiating (19) and substituting £ = 1 according to the second 
condition gives 

an = -X (IX) 

According to the third, condition, 

m = 0.14 (III) 

X X 

Therefore, a = -0.14 and n = , Since the value of X is 

a 0.14 


usually of the order of 0.43 to 0,44 (according to the resistance law 
derived for the pipe), the approximate result is: 
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D» 3 


and the required function 


B 


= f( 0 


will he a cubic parabola 




By substituting the obtained relation in the fundamental equation (18) 
there is obtained a differential equation the solution of which shot) Id 
give <p = cp (I), that is; the curve of velocity distribution: 


2 r _ 

— / cp 2 rd£ 

e 2 J 
* o 


_ G - + X ^ (1 - ^ 3 ) 2 I? 

2b 2 <W 0 I 


( 21 ) 


3. Velocity Distribution 


Integration of the obtained equation in finite form is not possible 
For this reason it is necessary to proceed to an approximate integration 
By integrating the inertia integral on the left in the form of a power 
relation and setting 


f; 



cp a sdt = J 


it is assumed that 


J 


% 


r k 

+ m x f 


(22) 


'The constants a x and m x will be determined from the boundary condi- 
tions for J. For 1-0, 


cp 2 c 

lim J = lim — - 

p — > 0 {: — 2 1 



2 


for f = 1 
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J 


l 



o 


cp 2 td£ = constant = J 0 


Then 


and 


a i 


& x + ma 






Elx 





£ 


k 


It la practically more convenient in the preceding expressions to pass 
from the universal velocities cp and <p m to the nondimens ional 


u cp 

^ c fk 


which are directly measured in the tests. On dividing the above expres- 
sion by cp 2 and letting 


I 




( 

\ 

\ 


u 

u m 


2 

t 

/ 


Id* 


there is obtained 


I 




t k 


Substituting in equation (21) the expression for I and solving for the 
first derivative of the nondimens ional velocity yields 
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*P _ 3 

dS = ~ X 




l\,k 

■ i j 8; 



(23) 1 


According to expression (ll) 

03 c 
2b" 


* 2-0 o <Pm *o - 1 


Substituting this expression in equation (23) yields 


dcp 

dF 



/ — — — — ’ — 
/*o*m (1-W 0 )(1 



) + 1 



(24) 


The expression ^ o %^ ( 1 - 2J 0 ) completely characterises the flow state 
in the diffuser, since it- takes account both of the geometric ('3 0 ) and 
the dynamic parameters cf^ and I Q depending on the Reynolds number. 

Let 


WU “ 2I 0 ) = D 


Kzr 


and denote this nondimensional magnitude as the "diffuser parameter." 

It is not. difficult to see that for the "zero diffuser" (pipe) D = 0. 
For the case of flow of an ideal fluid in a diffuser, that Is, when the 
velocity does not depend' on 5 and at each, section u = u m . 


l 



o. 


1 The minus sign before the square root is chosen from the consider- 
ation that over the entire range of variation of I the derivative , 

dcp 

— is less than zero. 
d| 


linn 


III ■■ I III! 


i i an 


ilia 
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In the case of a. diffuser and a converging pipe with turbulent flow 

where in computing I- r the function — entering the integral ie al— 

% 

ways less than unity, the result is 


I o< 


1 

2 


for which reason D > 0 for the diffuser and D < 0 for the convergent 
pipe (^ D < 0). Integrating (24) with respect to I gives 





/\ S D(1 -S k ) + 1 
l-l 3 


d| 


(25) 


Nothing, as yet, has been said regarding the value of the exponent 
k in the integration of the inertia integral. It is convenient to pro- 
ceed as follows: Assuming any value of k, to compute the integral 

(25); then, having the relation — = f(t), to' set up the values 

% 

I = f(£) and choose a new value for k in better agreement with the 
foregoing relation. By repeating this process several times it is not 
difficult to arrive at a value of k which best satisfies both expres- 
sion (25) and the approximation I = f(£). By this procedure, however, 
it is necessary to compute the integral (25) graphically, since finite 
computation ie possible only for a few values of k. 

In attempting to obtain the velocity distribution formula in finite 

3 

form, it was assumed that k = — which, as will be 3een later, is in 

2 3 

good agreement with all the conditions discussed. For k = ^ the in- 
tegral (25) is obtained in finite form by the substitution 



x 


leading to the integral of a rational fraction. Without going into de- 
tail, this gives the final result: 
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<Pm - ? 


1 / lQa (x + !)(/ gTI - 1) 
* ' (x - 1) (/ D + 1 + 1) 


2 ,/2D + 1 


tanh 


x 

/~2D + 1 


—1 

— tanh 


D + 1 
2D + 1 



(26) 


The proposed method of approximating the mixing length function by 
a cubical parabola gives a very simple law of velocity distribution for 
the straight pipe. Setting in expression (25) B*0 and carrying out 
■the integration gives for the pipe 


<Pm 


- cp * ~ tanh 1 

A 


( 



(27) 


.This very simple formula, as will be shown in part II, is excellently 
verified by experiment. 

It is of interest that the absolute character (independent of the 
Reynolds number) of the velocity distribution profile in the form 

<P m - <P * f(t) 

in the case of the pipe doeB not hold for the diffuser, for which case 
the "diffuser parameter" depending on the Reynolds number enters the 
velocity distribution (formula (26)), 

4. Resistance Law 

To make use of the foregoing derived formula for the velocity dis- 
tribution in the diffuser, it is necessary to be given the diffuser 
parameter 

D - VPm 2 (l ~ 2I o) 

that is, the values I 0 and <p m . The value of the inertia integral 
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yill be known, since the velocity profile is known. It. is suggested 
that an approximate value for Io first be assumed and then, having 
obtained the profile (26), to correct this value. Generally the succes- 
sive approximations converge very rapidly, since the values of I Q 
depend little- on the shape of the velocity profile. 

% 

It remains, for computing D, to assume the value m_ « — , which 

* 

must be associated with some characteristic of the flow, for example, 
the Beynolds number. ThiB means that it is necessary to find the re- 
sistance law. 

To obtain the resistance law, the solution obtained must be con- 
nected with the laminar layer near the wall. Since in the laminar layer 
only the effect of the viscosity should be assumed, consider where these 
equations lead to if the viscosity friction is considered instead of the 
turbulent friction. 

By retaining for the present the assumption of, radial flow (v n 0) 
also for the case of the laminar sublayer it is found that the intensity 
of the laminar friction is expressed as 


lam “ r 


du 1__ dt 

dfl r 3 as" 


(28) 


The assumption of radial flow leads to the result that the laminar fric- 
tion, in contrast to the turbulent, varies inversely proportionally to 
the third, and not the fourth, power of r. 


The equations of the laminar motion on the assumption of radial 
flow, on the basis of the fundamental equations and the above expression 
for the friction, become 


2t 


2 


5 

r 


l 8p J_ [" aft _ i at " 

p c)r r 4 -&a S 8 afl _ 


1 dp 2V 
p Sd ' " r 3 


dt 

08 
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On eliminating the pressure p "by the method of cross differentiation 


there is obtained after dividing by 



4 t dt v d 3 t s dt d 
r d-0 ~ d# 



-In the foregoing expression the right side does not depend on r; whereas 
the left side is inversely proportional to r. This absurdity clearly 
shows that the viscous fluid cannot flow radially in a conical diffuser. 
With this result. It becomes evident that radial flow in a turbulent 
fluid is, strictly speaking, impossible because the internal friction in 
a turbulent flow is the sum of the turbulent and the viscous friction. 

It is not difficult to see, however, that the deviation from radiality 
in the main body of the flow should be negligible, being of the same 
order of smallness as the proportion of the viscous friction is to the 
total friction. In other words, the flow in the diffuser may be thought 
of as a strictly radial, "purely turbulent " flow with superposed small 
nonradial disturbances arising f^'om the viscosity. 


The greatest deviation from radiality should be expected near the 
wall where the proportion of the viscous friction is particularly large. 

In order to form a very approximate picture of the flow in the laminar 
sublayer the following reasoning is given. It has been seen that for 
the radial flow assumption valid in the main bod,y of the flow the total 
friction should, on the basis of the Reynolds equations, vary inversely 

4 

proportionally to r . On the other hand, on the basis of the Prandtl 
formula and the assumption of universality of the none? I mens! onal mixing 
length curve the purely turbulent part of the friction should likewise 

4 

be inversely proportional to r . This justifies the assumption that 
the laminar part of the friction should have the same property with the 
same degree of accuracy. With the object of obtaining a smooth profile 
of the friction distribution, it is permissible to extend this property 
also to the purely laminar layer, that is, assume that 


T lam 




1 

~r 

r 


In this way, the assumption ma.de in solving the problem of the turbulent 
flow in a straight pipe is extended. It was assumed there that at the 
wall there is a viscous layer the motion in which is subject to the 
equations of a viscous fluid, but the value of the friction fib the wall, 
entering as a boundary condition, is determined by the turbulence law 
(resistance formula) . Here, by extending the foregoing assumption, 
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it is Considered that not only the value of the friction at the wall hut 
also the law of the distribution of the friction along the diffuser is 
determined by the main turbulent part of the flow. 

In the laminar layer, the intensity of the frictional stress in 
terms of the derivatives of the 1 velocity in spherical coordinates is 
expressed: 


T laiu ~ 




On substituting, on the basis of the above assumption, T ^ ara = , 

it is found that for the obtained expression to be independent of r, 
it is necessary that in the laminar layer the velocity components u 
and v be inversely proportional to r 3 , a condition which also does 
not contradict the continuity equation. Thus, near the wall there 
should occur a somewhat more rapid decrease in velocity along the dif- 
fuser than in the regions with rapid flow. With the velocity near the 
wall decreasing more rapidly, the direction of the component v, from 
considerations of continuity, should be from the wall to the axis; that 
is, the direction of a line of flow near the wall deviates from the ra- 
dial direction inward. 

The angular deviation of the lines of flow from the radial direc- 
tions (very small, of course, in absolute magnitude) near the wall 
where the effect of the viscous friction is large may be comparable 
with the small divergence angle of the diffuser. In this case it is to 
be expected that this deviation may affect the main turbulence of the 
flow in which the mean velocity vectors deviate somewhat from the ge- 
ometric radii. The problem of the deviations from the radial direction 
in part III of the present paper will be considered later. 

It may be well, now, to proceed with the direct derivation of the 
resistance formulas by the general method of considering the velocities 
at the boundary of the laminar layer. The velocity distribution 
u * f(-9) in the laminar layer is assumed to be linear: 


u 


( 


Su ^ 
Sy 


^y-o 


y 


where y is the distance from the wall. This formula corresponds to 
keeping the first significant term in the development of the function 
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u = f(y) 

in a series, a procedure which is entirely permissible in view of the 
extremely small thickness of the layer. Substituting 


y » r(d 0 - «) = r0 o (l - |) 


and expressing 



through the frictional shear at the wall: 


( 





gives 


u 


h 1 2 * * * * -8. 


r 8 v 


(1 - 0 


(29) 


The thickness of the laminar layer is assumed to follow the well— known 
law of Karman: 


6 


a 


v 


v 


* 


where, according to Karman’s assumption, a is a universal constant. 
The value ? = |* corresponding to the edge of the layer will then be 


- 1 - 


rfl. 


- 1 


a 


Vr 

bC 


1 " 

As K. K. Fediaevsky has shown; the assumption of the constancy of 

a corresponds to the assumption of the constancy of the critical 

Reynolds number computed for the thickness of the layer. There is 

therefore every reason to expect that in passing from the pipe to the 

diffuser the constant a does not appreciably change, 
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The nondimens ional velocity 
talned from equation (29) • 



at the edge of the layer is oh— 


9 5 = 




(1 - ?*) 


a 


The equations 


and 


<P ~ *^8 = a 


I 


_ ' s * 



will he substituted in the formula for the velocity distribution (26) . 
First of all, to find x - x* corresponding to ? = S *, apply the- ap- 
proximate expression for the square root as is always possible because 


of the extreme smallness of the magnitude 
There is obtained: 



compared with unity. 



By proceeding similarly with the power of I* there is obtained 


and 


•= 3 aVr 

V' 3 » 1 - - 

2 


+ i Dow 


4 R 


In substituting x = x* in expression (26), the term — — in com- 

4bd 0 - 

parison with unity will be neglected .throughout except, of course, 
term x — 1 . There is obtained: 


the 
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_ 1|, u "o , 3Da(/ X) + 1 + 1) 

«Pm - a * - I log log L 

X L Vr 8(/ D +1 - 1) 


f -rl 

— 2 »/ 2D + 1 ^ tanh 


1 -tax *- 1 V 


/ 2D + 1 


X 2D +1 y J 


The resistance coefficient and the Reynolds number are introduced as was 
done in the solution of the problem of the flow in a pipe . In analogy 
with the case of the pipe, let 


and 



= 


u m^ 

V 


where R = r9 is the local radius of the circular cross-sections of 
o 

the diffuser. This gives 




% 


P% = j/2 

T o /T t 


ba ^ v ,.r 2 d 
o _ ' r o 

v r v r 



B m/ C f 

/1> 


By substituting in the preceding expression and' passing from hyperbolic 
functions to logarithms the resistance law for the diffuser is obtained: 
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i°e (Bb /o f ) + a 
/% X /fi ■ 


where 



1 


+ / 2D + 1 log 


(l + / 2D + 1)(/2D +1 - /D + 1) 

(/ 2D V i - i)(/sd' + i + /jj'+'i) 



(30) 


On passing from the diffueer to the straight pipe, that is, setting 
D = 0, there is obtained on evaluating the indeterminate expression by 
the rule of L* Hospital: 


where 



—i— log (1^ /C f ) + A 
V. /? 


A 



log 

X /!> 


4 


(3D 


As may he seen, the form of the resistance expression for the diffuser 
and that for the straight pipe is identical. It is very interesting 
that the diffuser parameter enters only into the makeup of the free 
term A and has no effect whatever on the terras containing the resist- 
ance coefficient and the Reynolds number. This is not entirely 
unexpected, because the resistance formula results from the shape of 
the velocity distribution curve near the outer boundary of the laminar 
layer, where the effect of the inertia forces is negligible and the 
mixing length, as in the case of the straight pipe, increases propor- 
tionally to the distance from the wall. 

The numerical value of the free term A usually is determined 
experimentally for the pipe. As is known, the value of A is abso- 
lute, that is, independent of the Reynolds number, and is of the order 
of 4 to 5. 1 


^According to the tests of Itfikuradse, A = 4.88: according to the 

present tests with the straight pipe A = 4.01. 
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There is every reason to expect (since the form of the expression for 
the resistance is the same as for the pipe) that the absoluteness of A 
extends also to the diffuser, while a variation is possible in the 
magnitude of a which for different values of the parameter D will 
be obtained as a solution of the transcendental equation 

A = f(a, D) * constant 

The circumstance that the diffuser parameter D, which takes account 
of the effect of the inertia forces, does not enter the variable terms, 
of the resistance formula leads to the conclusion that the form of the 
expression is not due to the approximation used above for the inertia 
integral. If it were possible to solve accurately the fundamental dif- 
ferential equation, the same form of expression would be obtained for 
the resistance formula. 

The approximations assumed for the mixing length and inertia in- 
tegral show up only on the free term. For this reason, no far reaching 
conclusions will be made with regard to the dependence of a on D 1 
or physical explanation of this dependence, particularly since for 
practical computations not a but A is required., which in any case 
is obtained from experiment. 

It is necessary to make one more remark with regard to the resist- 
ance formula. The Reynolds number entering the formula 

pje = ^ 1Jl m _ r ^o _ '^o 


is variable along the diffuser length. The 
ever, according to the assumption of radial 
the diffuser length: 


C f 



2fp r 4 
4 2 

r pt^ 


friction coefficient, how- 
flow should not vary over 

= gfo 

2 

Ptm 


This lack of correspondence is obtained as a result of the foregoing 
assumptions in considering the viscous friction at the wall which dis- 
turbs the radial character of the flow. The error thus obtained is 
vanishingly small. For small divergence angles of the diffuser the 
Reynolds number Re m varies so little over the diffuser that, Cf, 


*It is seen, for example, since D = f(Re) that also a in the 
diffuser will depend on the Reynolds number, that is, no longer be an 
absolute constant. 
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depending on the logarithm of Be m , undergoes quite, a negligible change. 


II . EXPERIMENT AND COMPARISON WITH THEORY 


In order to check all the assumptions made and results obtained in 
the preceding part, a series of tests were carried out on two conical 
diffusers with divergence angles (angle botween axis and generator of 
cone) of 1° and 2°. In comparing the experimental with the theoretical 
results it was assumed that the comparison of the velocity and pressure 
distributions was of greatest importance. For this reason the experi- 
mental determination of stich secondary factors as the frictional dis- 
tribution over the cross section of the diffuser, the. distribution of 
the mixing length, and so forth, were entirely omitted on the assump- 
tion that agreement or disagreement between theory and experiment for 
the main factors implies also the corresponding condition for the 
secondary factors. 


(a) TeBt Setup 

The diffusers on which the tests were conducted were of plywood 
construction having sheets, curved into conical segments, which'were 
attached by wooden rings the Internal diameters of which as far as 
possible were determined by accurate computation. The wooden rings 
were connected to each other by stringers disposed over the generators 
of the cone. The entry opening of the cone had a diameter of 2 ho mil- 
limeters and was connected by means of a transition piece with the 
cylindrical pipe previously investigated (reference 1) having a length 
of about ?0 calibers. The diffusers were each 6 meters long. Behind 
the diffusers was placed a 2-meter section of a straight pipe having 
the diameter of the outlet section. This section was followed by a 
short accurate converging pipe that reconducted the flow into the pipe 
of 2 4C— millimeter diameter. The latter pipe with the aid of two elbows 
was connected to the suction chamber. The inner surface of the diffuser 
and of the transition lengths was covered with shellac . The entire ap- 
paratus was suspended from the ceiling of a long passage of the lab- 
oratory. The lines in the diffuser cross sections along which the mean 
velocity distribution was measured lay in the hori zontal plane passing 
through the axis of the apparatus. The first section was located in 
the center of the forward transition piece. The distance x, in 
meters, between the succeeding sections measured from the fihst is 
given in table I . 
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TABLE I 


Section 

B 

2 

3 

| 

5 

1 

6 

D 

8 

9 

. 10 

11 

■a 

13 

X 

0 

H 

2 

2.2 

1 

2.6 

2.8 

B 

3.2 

fl 

3.6 

3.8 



In distributing the sections along the diffuser it vas assumed that 
over a certain portion beyond the transition pipe no typical diffuser 
flow will be observed to which the assumptions made in part I were ap- 
plicable. Over this portion there will be a gradual transition from the 
conditions of flow in the straight pipe to the conditions of flow in the 
diffuser; that is, there will be a kind of diffuser entry length 
(Anlaufstrecke) . On the basis of this fact, the main sections 3 to 13 
were placed in the middle part of the diffuser. The test, as will be 
shown, justified these assumptions. 

The velocity distributions were measured in sections 3 to 13 with 
the aid of the same pitot tub6 used In investigating the velocity dis- 
tribution in the cylindrical pipe and the same coordinate apparatus 
(reference 1) . Since the openings, through which the stem of the tube 
passed into the diffuser and came out through the opposite side, were 
made along the normal to the inner surface of the diffuser, the tube 
stem was slightly curved. It thus seems reasonable to assume that the 
velocity distribution was measured over the arcs of circles of greater 
radius, approximately corresponding to the arcs measuring the angle d. 

To measure the static pressure distribution along the diffuser, the 
usual system of brass tubes was mounted along the bottom stringer. The 
tubes were placed 200 millimeters apart from section 1 almost to the end 
of the diffuser. 


(b) Tests 

The velocity distribution in the chosen sections of the diffuser 
was measured for a single maximum discharge rate of air. The velocities 
were measured in each section from wall to wall, in order to be able to 
judge to some extent the axial symmetry of the velocity profiles. The 
points at which measurements of the velocity were made in each cross 
section were more often near the walls than near the diffuser axis . In 
testing, two micromanometers were read, one connected with the pitot 
tube (h^), and the other a control (h c ) connected with a static 

tube placed in the cylindrical pipe ahead of the diffuser. After trav- 
ersing each section of the velocity profile, the pitot tube was placed 
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on the diffuser axis and the air discharge was varied from zero to the 
maximum with the aid of the iris diaphragm at the chamber.- Simultaneous 
readings were taken on the velocity (h^) and control (h Q ) manom- 
eters, and the magnitude was computed: 



(32) 


which for each section was plotted as a function of h c . It was thuB 

possible to compute the absolute value of the velocity at the diffuser 
axis as a function of the reading on the control manometer: 


Uj£ = / — F v Ky A h c 

/ P 


where 


P 


0.125 


g88 Patm 
760(273 + t) 


is the density of the air with corrections on the temperature of the air 
and the atmospheric pressure; 


F v and k^. sine of angle of inclination of manometer tube and 
coefficient of manometer, respectively 

5 T coefficient of pitot tube 


7 


density of spirit in manometer 


Thus, by instantanecrus reading of two manometers, it was possible 
to determine the velocity at -any point of any section as a fraction of 
the velocity at the diffuser axis in the same section: 


u 


u. 


m 


/£ 

/ Ahc 
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In determining the velocity distribution by this formula, near the 
diffuser walls a correction was made for the effect of the walls on the 
readings of the pitot tube. It is known that the presence of a wall 
near the stem of the tube gives rise to a certain increase in the ve- 
locity between the stem of the tube and the wall, and this leads to a 
drop in the static pressure and results in an increase in the manometer 
readings above the true values. The correction was applied by multi- 
plying the radical on the right side of this formula by the correction 
coefficient £ . The dependence of this coefficient on the distance y 
from the wall to the center of the front opening of the tube is shown 
on figure 3. This curve was obtained from tests in the straight pipe 
(reference 1) where the velocity distribution near the wall was measured 
both by the pitot tube and the total pressure micro nozzle, on the read- 
ings of yhich the effect of the wall was vanishingly small. The test 
points on figure 3 were obtained for various Eeynolde numbers. As may 
be seen, with increasing distance from the wall the values of £ rap- 
idly approach unity. 

In addition to the correction near the wall, a correction was made 
for the change in diffuser cross section caused by the stem of the tube 
in determining the absolute value of the velocity u m on the diffuser 

axis. This correction coefficient u was determined as: 



2 2 

where S = jtr $ 0 circular area of the diffuser cross section, 

Sjjjp = 2rfi 0 d part of diffuser cross section taken up by stem of tube 

(d diameter of tube). The justification for applying the correction 
by this method is given in a paper on tests on a pipe (reference l) . 


The distribution of the velocity ratio — for the various sec— 

u m 

tions of the diffusers of 1° and 2 ° is shown in figures 4 and 5, where 
the magnitude 


E =• 



z 



is laid off on the axis of abscissas. The values z and z 0 , the 

lengths of the arcs of the circles having their center at the vertex of 
the cone and passing through each section, were read off directly with 
the coordinate apparatus by divisions carried on the stem of the pitot 
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tube. On these figures the velocity curves are shifted upward with re- 
spect to each other. The dots and circles denote the measurements at 
the right and left halves of the velocity profile. 

In examining these curves there is observed first of all that: (1) 
The symmetry of the profile with respect to a vertical plane is suffi- 
ciently good (the dots and circles are not very far removed from the 
mean curve) , (2) the assumption in part I on the similarity of the ve- 
locity profiles for sections 3 to 13 is observed to be well satisfied. 
An exception occurs in the case of the 1° diffuser for section 5 : which, 
for an unexpected reason, drops out of the general series. 

Of interest on the above curves are the velocity distributions for 
sections 1 and 2. The latter distributions (especially for the 2° dif- 
fuser) clearly indicate the transition region from the straight pipe to 
the diffuser. Eere, it should be remarked that the velocity distribu- 
tion in section 1 is not characteristic for the pipe and is also 
transitional, since it must be supposed that the effect of the diffiiser 
extends somewhat ahead into the pipe. On the basis of the assumption 
that the shape of the velocity profile is entirely determined by the 
static pressure gradient, it may be said that the region of influence 
of the diffuser on the pipe will extend into the latter up to a point 
where the pressure gradient in magnitude and sign will agree with the 
corresponding value for the pipe. 

After measuring the velocity profiles over the diffuser sections, 
pressure distribution measurements were carried out with the aid of the 
above-mentioned system of static tubes. The measurements were made lor 
five discharge rates for the 1° diffuser and eight discharge rates for 
the 2° diffuser. 

The results are given on figures 6 and 7, where the lowering in 
pressure in kilograms per square meter is plotted as a function of the 
distance from section 1. As may be seen, in both diffusers there was a 
positive pressure gradient. The static pressure distribution trans- 
verse to the diffuser was also measured. Ho pressure gradient with re- 
spect to ■d was observed. This completed the test program. 


(c) Checking of Fundamental Assumptions 

The confirmation of the fundamental assumption made in constructing 
the theory, namely, the assumption of radial flow will be discussed 
first. On applying this assumption to the equation of continuity, the 
result (formula 5) that the product ur 2 should 1 not be a function of r 
was obtained. The values of this product on the diffuser axis where 
u = u m will be computed. In tables 2 and 3 are given the values of 
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and r for both diffusers for maximum discharge rate, and the products 
are -computed. The values of u m for the various sections of each dif- 
fuser, with the aid of the curve A = f(h c ), was reduced to a single 

reading of the control manometer, that is, to a single discharge rate. 


TABLF II .- VALUES OF r, AND u m r 2 FOR 1° DIFFUSER 


Sec- 

tion 

3 

D 

5 

6 

n 

8 

9 

... ... 
10 

11 

12 

g 

B 

15-73 

15.93 

16.13 

16.33 



16.93 

17.13 

17.33 

n 

17.73 

(m/ sec) 

21.0 

20.2 

19 -9 

19.75 

19.1 

1.8.7 

18.5 

18.05 

^ b 
c— 

r-i 

17.34 

16.9 

V 2 

5200 

5130 

5190 

5270 

5225 

5225 

5310 

5300 

5320 

5340 

5310 


TABLE III.- VALUES OF r, u^ AND U n r 2 FOE 2° DIFFUSES 



0 

J 

■ 

5 

6 i 7 

\ 

8 

9 

10 

— 

11 

— 

12 

13 

r 

(m) 

. 

8.88 

9.08 

9.28 

! 

9.48. 9.63 
1 

9. 88 

10.08 

10.28 

10.48 

10.68 

10.88 

% 

(m/sec) 

18.30 

0 

CO 

« 

H 

17.18 

1 

16.65 |l6. 18 

1 

15.50 

14.95 

14.50 

13.66 

13.17 



12.63 

2 

1443 

1470 

1480 

IBM 

B BB B 


1525 

1540 

1500 


1500 


For greater clearness the values of the product i^r 2 are plotted 
aB a function of x in figures 8 and 9« As may he seen, the product 
Uj^r 2 does not remain constant hut increases in a regular manner along 
the axes of the diffusers. The regular character of the increase in 
Uj^r 2 indicates that the obtained result is not due to experimental er- 
ror hut to a certain regular disturbance in the radial character of the 
flow that arises from the effect of the viscosity. The difference 
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■between the maximum and minimum values of u m r 2 is, however, sufficient- 
ly small and does not exceed k percent for the 1° diffuser and 7 percent 
for the 2° diffuser. The regularity of the increase in u^r 2 is broken 
only in the last sections of the 2° diffuser, where evidently some other 
unaccounted-for effect enters. It is possible that the irregularity at 
these points may be due to some error in the test setup. 

Thus the radiality of the flow at the diffuser axis (and therefore 
also over the entire region of the flow) is satisfied only with a cer- 
tain, though rather large, degree of approximation. A more detailed 
discussion on the radiality of the flow will be given in part III of 
this paper. For the present, however, the obtained increasing values 

u m r = f(r) will be replaced by me an^ constant values (u m r ) m . On 
figures 10 and 11 the values of { u m r ) m are plotted and the mean 
values indicated for various discharge rates for both diffusers. 

Another method of checking the radiality of the flow is to study 
the increase in the static pressure along the diffuser. In the first 
part the result was obtained that in the case of radial flow the gen- 
eralized pressure p - t should vary along the diffuser according to 
the law 

Pi * - ~ ~ + ( Pi Q ~ —j- ) (formula 8) 

4 r kr c '' 


At the wall, that is, for ? = 1, the pulsations should die down 
( T rr = 0) and Pi = p. Tiius, in the case of the correctness of the 
foregoing law of Increase in pressure, there should in plotting p as 

a function of ~ be obtained straight lines, the slope of which is 
r 

determined by the value of G: 


k dp 
p d(r 4 ) 


(33) 


This plot for the two diffusers is shown on figures 12 and 13. As may 
be seen, the law expressed by formula (8) is excellently confirmed. 

The various straight lines on these figures correspond to the values 
(u m r ) m on the previous figures. For each straight line is indicated 
the value of G determined by formula (33). 
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In measuring the pres&ure distribution for ,the 2° diffuser use was 
made of a more improved system of manometers .as compared with that used 
for measuring the pressure in the 1 ° diffuser. Because of this it was 
possible in figure 13 to plot the test points directly, while in plot- 
ting figure 12 the values of the pressure had to be taken from the 
comparison curves (fig. 6). This explains the greater scatter of the 
points on figure 13 as compared with figure 12. 

On these curves the values of x corresponding to the sections at 
which the velocity distribution was determined are indicated below. It 
is of interest that the pressure variation law is confirmed over a con- 
siderably greater length than the law of similarity of velocity 
profiles. This should also be the case since according to equation 
(11) the value of the magnitude G characterizing the change in pres- 
sure depends not directly on the velocity distribution but on a certain 
integral of this distribution which should not be -very sensitive to a 
change in the velocity profile. 


(d) Check of the Velocity and Resistance Laws 


On figures 14 and 15 are given the comparison of the experimental 
and theoretical profiles ~ = f(?) of the velocities determined by 

Um 

formula (26) for the maximum Reynolds numbers attainable.’ 1 The value 
of the diffuser parameter D for the computations by formula ( 2 6) was 
determined to a first approximation by the formula 


D = * 0 <P m (l ~ 2 I 0 ) 


where the value of I Q was computed with the aid of graphical integra- 
tion of the mean experimental velocity profile . After obtaining the 
theoretical profiles shown on figure 14 and 15 , the value of D was 
corrected Hy- determining I Q from the computed theoretical profiles. 

The second approximation for I> differed so slightly from the' first 
that further approximations of the velocity profile seemed superfluous. 
The value of entering the formula' for I) was determined as 


cp m ■* 




/ 





Rote the transition here from the value s cp^ — cd to the values 


— , determining cp according to the resistance law. . 
u m 
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from the derived resistance law ( 30 ) where the value of the absolute 
constants X -and A, on the basis of teats with cylindrical pipes, 
were taken as 

X = 0.434 

and r 

A = 4.01 

The Reynolds number Re m was determined as the mean of the 
Reynolds numbers 


R©m 


r ’ 3 o ll m ‘ 3 o^ u m r " 
v rv 


over the diffuser sections. The possibility of this averaging was jus- 
tified by the fact that the change in Re m along the diffuser 
practically had no effect on the values of and cp^. 

The final values of D and for both diffusers for mean 

Reynolds numbers were obtained as: 

For the 1° diffuser Re m av = 178,000 cp^ = 26.65 R = 2.27 

For the 2° diffuser Re m av = 177,000 = 26.55 R = 5.5 6 

According to figures 14 and 15 it may be stated that the theoretical 
formula satisfactorily agrees with the test results. Some deviation is 
observed in the middle portion of the curve where the test points have 
a tendency to drop below somewhat, The reason for this deviation must 
be sought, apparently in the fact that the actual distribution of the 
mixing length in the diffuser differs somewhat from the cubical parab- 
ola assumed. 

To illustrate the degree of accuracy of the solution of the funda- 
mental differential equation by approximating the inertia integral by 
the exponential relation 


3 



figures l6 and 17 are presented, where the continuous curve gives the 
above integral and the points give the results of computing I by 
means of the graphical integration of the obtained theoretical curve of 
the velocity distribution. The discrepancies between these values of 
X are so small that it may be affirmed that the degree of accuracy of 
the solution of the fundamental differential equation by this approxi- 
mation is sufficiently high. 
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the velocity distribution. The discrepancies between these values of 
I sire so small that it may be affirmed that the degree of accuracy of 
the solution of the fundamental differential equation by this approxi- 
mation is sufficiently high. 

It is very unfortunate that formula (26) does not enable the an- 
alytical computation of the values of the inertia integral, since the 
problem reduces to incomputable quadratures. Otherwise by expressing 
I Q in terms of D by means of the above integration the dependence 

of D could be found only on the Reynolds number and the divergence 
angle of the diffuser. 

As an illustration of formula (27) of the velocity distribution 
for the straight pipe (zero diffuser), figure 18 shows the results of 
the tests in the pipe as given (reference 1) and, as may be seen, for- 
mula (27), for the values of X chosen on the basis of the resistance 
law, is in excellent agreement with the test results. It is in- 
teresting to note that the curve representing formula (2?) lies very 
near the curve drawn according to the theory of Mattioli, which until 
now of all semi empirical theories best agrees with experiment. This 
consideration is of great importance since formula (27) is extremely 
simple; whereas the theory of Mattioli does not lead to any finite for- 
mula at all, and to draw the velocity distribution it is necessary to 
carry out a graphical integration. 

For checking the resistance formula in the usual manner by com- 
paring the theoretical and experimental laws of variation of the 
resistance coefficient with the Reynolds number, it would have been 
necessary for to find the experimental values of the frictional shear 
at the diffuser wall. It was considered illogical to carry out the 
above operation, since It would then be necessary again to make use of 
the theoretical equation of motion that gives the relation between the 
friction, the inertia forces, and the pressure gradient. For this rear- 
son it waB considered of greater interest to compare the theoretical 
value of the magnitude G, characterizing the increase in static pres- 
sure along the diffuser, with the corresponding test values obtained ty 
measuring the slopes of the curves on figures 12 and 13 . 

The theoretical value of G may be obtained from equation (ll): 


G 



29, 


r 

1 


<p 2 |di 


-xV 


2b a cp 2 


m 



_JL 





1 


\ 

"\ 


) 


(34) 
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where is no other them the value (u m r ) m , depending on the rate 

of discharge of the fluid through the diffuser, that is, on the 
Reynolds number 



\ 




By integrating the theoretical velocity profile, I Q is determined; 
cp^ is determined as a function of the Reynolds number according to the 
resistance formula. 

Thus, expression (34) is the theoretical dependence of G on the 
Reynolds number while expression (33), which determines G by the 
slopes of the test curves, is the corresponding experimental dependence . 

On figures 19 and 20 are given the curves of the theoretical and 
experimental relations G = f(t m ). As may be seen, in both cases the 

test points lie below the theoretical curves . Part III of the present 
paper will show more in detail the reason for the obtained (not very- 
large, it is true) discrepancies in the values of G. 


In concluding this part, the theoretical curves of the friction 
distribution transverse to and along the diffuser are given. By formula 
( 10 ): 


■S G \ / d G \ 

" 5 ( 23 0 J - j 

Oh 2 / \ Oh 2 


By substitution, according to expression (ll), 


■8 G 

- EVo - 1 

2b 


there is obtained 

— “6 [2^8 0 (J - J D ) + 1] 

TO 

According to the assumed approximation for the inertia integral. 
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T 




- 2J„ 


) 



g 

By substituting J Q *= cp And introducing the diffuser 'parameter D 
there 'is obtained finally 


T 



= £ 



(35) 


On figure 21 are constructed the curves — = f(6), according to 

T o 

formula (35) , for the diffusers and cylindrical pipe. It is interesting 
how such small diffuser divergences as 1° and 2° , which seem insig- 
nificant at first sight ; affect the character of the friction distribu- 
tion. As may be seen, near £ =? 0.6 the local frictional shear for 
the 1° diffuser 1 b approximately twice and for the 2° diffuser about 
four times the value for the straight pipe . 

On figures 22 and 23 are drawn the friction distributions along 
the walls of the diffusers according to the law 



4 

r 


for various values of = (u^r ) m . The values of f Q is expressed 
as 


f o - ~ ^Pt 2 
O g 1 m 


t \2 

) 

9m. ^ 


III . SOME CONSIDERATIONS WITH REGARD TO TEE DEVIATIONS 
FROM THE RADIAL FLOW ASSUMPTION 


It has already been pointed out that the none one tancy of the values’ 
of u m r c along the diffuser axis indicates a systematic deviation from 
the assumed radial flow. It is true that the deviation of the flow 
from the radial direction is extremely small. It is interesting. 
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however, to estimate to a first approximation what these deviations are 
and their effect on certain final results. 

Investigation of the deviation from the radiality is possible by 
considering the so-called ’'secondary flows," that is, the transverse 
components v of the mean velocity directed along; the circular arcs 
measuring the angles d. Hie direct measurement of these velocities 
is, of course, impossible. 

An attempt was made to use the continuity equation connecting the 
velocities v with the velocities u by analyzing with its aid the 
test curves of the velocity distribution. This attempt did not, how- 
ever, lead to any good, results, for it was found that the accuracy of 
the test for this purpose was insufficient . Some conclusions may, how- 
ever, be drawn from the obtained experimental lav of the variation of 
the magnitude u~r along the axis. 

It is assumed that at any point M with coordinates -8 and r 
within the diffuser there are two component velocities u and v 
(fig. 24). By adding u. and v vectorially there is obtained the 
modulus of the resultant velocity 

/ ~~2 2 

w = / U + V 


By prolonging the direction- of w to its intersection with the diffuser 
axis at point W, the point N may be approximately considered as a 
certain fictitious source producing a radial velocity field near point 
M. The distance of this fictitious source is denoted from the vertex 
of the cone by OH = Ar. By proceeding in the same way for each point 
within the diffuser, various values will naturally be obtained for Ar. 
In the general case, 

Ar = f(r, -3, u, v) 

The extent by which Ar differs from zero will characterize the degree 
of nonradiality of the true motion. 

This function is expressed through the coordinates of the point 
and the values of the velocities. From triangle OHM there is ob- 
tained (fig, 24) 

v 

= -d + tan — 
u 
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On account of the smallness of v it may he assumed that 
Bearing in mind the smallness of the angles, write 


tan 


-x v 




u 


u 


(r — Ar)-8i = rd 


whence, comhining with the ahov9 equation gives 


Ar = 


v 



v 

u + - 
d 


(36) 


The limit which this function approaches will he found as -3 approaches 
zero, that is, the value of this function on the axis of the diffuser: 


rum (j) 

•o —5* o 

Ar™ = lira Ar ; r- 

in. / n 

d — * O / v \ 

u m + lim [ - 

m _ o \ / 


From considerations of symmetry v = 0 on the diffuser axis. The limit 


of the ratio — will he found hy making use of the equation of con- 
’d 

tinuity, which for the two velocity components in spherical coordinates, 
will have the form: 


ou + 1 ov Su . v 

3r r od r rd 


Comhining the first and third terms of the equation results in 


5(ur g ) 

5r 


bv 

+ r — - + r 
dd 


V 

d 


0 


whence the required ratio 
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f 1 S(ur 3 ) Sv "N 




Sr od 


;) 


On passing to the limit there is obtained: 


lim ( 
o V 


S(%r ) 
Sr 


- u- ) 

4 o V Sd J 


(37) 


On the other hand., from L'Hospital's rule there is obtained: 

^ / Sv 


lim l — ) = lim . 
d — *• o \ -0 ' ■q -*-o V ^ 


(37a) 


Substituting this in expression (37) and solving the obtained expres— 

v 

sion for the limit of — gives: 

■a 


lim 
•a — » o 


/ v 1 

\ d ' 2r 


2 

°(%F_) 

dr 


Thus, 


hvT™ = 


Sfu^r 2 ) 

Sr 



S (u m r") 

Sr 


(38) 
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The value of the derivative ie experimentally determined from 

2 

the curves of figures 8 and 9 • The curve u^r « f(r) shown on these 
figures may, with a very large degree of accuracy, be represented as a 
linear function: 

u^ 2 « ar + b 


In determining the values of a and b, use was made of the 
method of least squares, which for the given case gives the generally 


formulas 

y» 

1: 

r3 

2u m r> 

n 

2' 

n 

Tl 

h — ■ 

2^ 

n 

2 r 

n 

2«»'-’ 

2“m r * 

n 

tl * 

2^ 

2' 

, U j 

2^ 

2' 


n 

n 


rt 

n 


2' 

n 


2' 

n 


n 



n 



where the sums are taken for n readings of u m * f (r) . (In the given 
case n is the number of diffuser sections in which was measured. ) 
For these cases there was obtained (for maximum Reynolds numbers): 

1° diffuser a « 88.6 square meters per second 

b * 3769 cubic asters per second 

2° diffuser a » 59.6 square meters per second 

b * 923*7 cubic maters per second 


This leads to the computation formula: 


Ar m = - 


(39) 


2 u„ — 


From expression (37a) It is found that the first derivative of the curve 
of transverse velocity distribution v with respect to d on the axis 
will be: 
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where the minus sign indicates that the velocities v are in the direc- 
tion opposite to those shown on figure 24, that is, from the wall to the 
axis. 


In tables IV and V are given the computed values of the function 
Zrfcj = f (r) according to formula (39). 


TABLE IV.- COMPUTATION OF Ar m FOB 1° DIFFUSER 


Sec- 

tion 

3 

n 

5 

6 

■ 

8 

9 


11 

r 

12 

13 

r 

(m) 

15.73 

15.93 

16.13 

16.33 

16.53 

16.73 

16.93 

17.13 

17.33 

17.53 

17.73 

% 

(m/sec] 

21.0 

20.2 

19.9 

19.75 

19. r 

18.7 

18.5 

18.05 

17.70 

17.34 

16.9 


-2.W 

-2.55 

-2.58 

-2.60 

-2.70 

-2.76 

-2.79 

-2.86 

-2.92 

-2.99 

-3.07 


TABLE V.- COMPUTATION OF hr m FOR 2° DIFFUSER 


Section 

.. 

3 

■B 

5 

6 

7 

8 

9 

10 

r 

(m) 

8.88 

9.08 

9.28 

9.48 

9.68 

9.88 

10.08 

10.28 

% 

(m/sec ) 

18.30 

17.80 

17.18 

16,65 

16.18 

15.50 

14.95 

14.50 


-2.00 

-2.06 

-2.13 

-2.20 

-2.28 

-2.40 

-2.49 

-2.57 


The minus sign before the values of Ar m indicates that the fictitious 
source is located not ahead of the vertex of the cone hut behind it. As 
may be seen, the values of £r m for the 1° diffuser constitute on the 
average 1 6 to 17 percent of the values of r and for the 2° diffuser 
23 to 25 percent of the corresponding values of r. This, as may also 
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be expected, is a consequence of the fact that the deviation from ra- 
dial! ty for the 2° diffuser iB greater than for the 1° diffuser. 


Moreover, it is seen that Ar m v aides over the diffuser length. 

The variations of &!&, it iB true, are not large in comparison with 
r and constitute over the entire range investigated 3-5 to k percent 
for the 1° diffuser and 5.5 to 6.5 percent for the 2° diffuser. 

The latter circumstance makes it possible to assume approximately that 
near the axis flow may actually be thought of as originating at the 
fictitious source behind the cone vertices at the mean distances 


and 


= -2.75 

for the 

1° 

diffuser 

2* a « -2.27 

for the 

£>0 

diffuser 


Since the true nonradial flow near the axis is approximately replaced 
by a fictitious radial flow originating at the second source, the equa- 
tion of continuity valid for radial flow may be applied to it; that 
is, write 



« 0 


whence 


u, 


m 



= = constant 


(**>) 


The correctness of this expression is readily proved by drawing curves 
similar to those of figures 8, 9, 10, and 11, the values u^r - Ar ra ) 2 
being laid off along the axis of ordinates. These curves are given on 
figures 25 and 26. On comparing with those of figures 10 and 11, it 
would seem that the constancy of determined by formula ( 40) is 

considerably better satisfied than without taking t he displacement of 
the source into account. In drawing these curves firm was considered 
as independent of the velocity on the diffuser axis, that is, independ- 
ent of the Reynolds number. This evidently corresponds with the facts, 
since in determining Ar m , formula ( 38) , both the numerator and the 
denominator may with a very large approximation be assumed proportional 
to the velocity 

For still greater assurance as regards the possibility of rep- 
resenting the nonradial flow as a radial flow originating from a 
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fictitious source, Ar m vas determined by still smother method. Let 
it be assumed that there is given a test distribution along the dif— 
fuser axis of the magnitude u^r . It may be asked what increment Ar m 
must be given to r in order that the vsilue of the product 
2 

i%n(r — Arm.) will deviate leaBt from any constant values for each 
Reynolds number. The statement of the problem- in this form leads to 
finding a minimum of the function: 




(r—Ar m y—u m (r — Ar m f 


where the first part represents the sum of the squares of the deviations 
of the values u m (r — Ar m ) from the mean taken for n readings. Re- 
moving parentheses, differentiating with respect to Ar m , and equating 
— to zero yields, after transforming the cubic equation in Ar m : 



The computations according to this equation, conducted with very great 
accuracy, gave for the diffusers for maximum discharge rate the roots: 
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1° diffuser 
2° diffuser 

The good agreement of the. obtained values of Ar m with those computed 

from the continuity equation very obviously indicates the acceptability 
of all the foregoing assumptions with regard to Ar. 

It is very interesting to inquire whether the results of the pres- 
ent investigation conducted for the region near the diffuser may be 
extended to the entire region of flow in the diffuser, assuming, of 
course, in the general case Ar = f(d). That Ar necessarily must 
depend on d follows from its very definition. Thus, for example, on 
the oasis of expression (36) at the wall, that is, for -3 = it 

it follows that: 


Ar m * —2.44 meters 
Ar m * —2.20 meters 


im Ar 


G — 1* 


zz - • 


o 

u 

V 


O V 
o 


U , 
“ + 1 


where 


lim 
v— » o 
u -* o 


u 


= lim 


u • 

v- 


3u 

3a 

3v 


° OG 


since on the basis of the equation of continuity 


3v f 

lim = lim j 
u-»o ^ u — o L 

V — » O V o 


1 3ur'~ v 

r 3r -6 


0 


while the value of 
Thus, 


3u 

3d 


at the wall is, of course, not equal to zero. 


(' t ' r )-g= J g n * 0 
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The function . Ar **. f (d) ~ evidently, .should be of the following form: 
It should increase sharply over the thickness of the laminar sublayer 
(where the effect of the viscosity, as has been seen, particularly dis- 
turbs the radiality of the flow) , and attain somewhere near the edge of 
the laminar sublayer -a maximum followed by a smooth drop in the direc- 
tion toward the diffuser’ axis to the foregoing numerical values. 

The possibility of representing the entire flow in the diffuser 
as radial with displaced source has, in the writer's opinion, a certain 
practical value. It is possible thus,' with a very high degree of ap- 
proximation, taking for the continuity equation the expression. 

b h»( r '~ to 



to' 'consider the phenomenon as subject to the equations of motion appli- 
cable to the radial flow but substitute in them throughout r — Ar in 
place of r. Thus there may be written instead of formula (6) 


and 


G 


1 °Pl ,5 

“ Tf TT ~ Ar) 

p d(r - Ar) 


( 6 * ) 


f = T(r - Ar) 4 


In place of expression (8) there is obtained: 


Pi 


Pi, 



1 

(r - Ar) 4 



( 8 ’) 


Expressions (33) and ( 3*0 determining the experimental and theoretical 
values of G become, respectively. 



P 



( 33 ’) 
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G = 2 



_ \ 2 "1 2 
) J 



(34*) 


These formulas were written for the axis of the diffuser where Ar = Ar m , 
still retaining the previously made assumption that 


&> = o and S£ * o 
3-b 


An attempt will he made to verify whether the agreement of the 
theoretical and experimental values of G is improved by the correc- 
tion for the nonradiality according to formulas (33*) and (34*). On 
figures 27 and 28, analogous to figures 12 and 13, are drawn the curves 






\ 2 

J J 


As may be seen, the introduction of the correction Ar m has almost no 
effect on the appearance of the above curves. 

On figures 29 and 30, a comparison is made between the theoretical 
and experimental values of G determined by formulas (33’) and (34'). 
As may be seen, comparing these curves with the previous figures 19 and 

20, it may be found that although the correction Ar m gave excellent 
results for the 1° diffuser, the test points for the 2° diffuser lie 
on the other side of the theoretical curve and at a greater distance 
than on figure 20. 

It would seem that it is not difficult to establish the reason for 
the result for bbe 2° diffuser. Up to now it has been assumed that the 
gradient of T rr along r is negligible in comparison with the cor- 
responding gradient of p, and this permitted consideration of the 
value of G as independent of $ . Evidently, although for the 1° dif- 
fuser this neglect was justifiable; for the 2° diffuser, where the 
intensity of the turbulence fluctuations should be higher, it is neces- 
sary to introduce a correction for T rr which is, of course, a func- 
tion of both ^ and r. The theory is as yet unable to contribute 
anything in this direction. 
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CONCLUSIONS 


The theoretical and test results presented on the problem of the 
turbulent motion in conical diffusers, are, of course, to a large ex- 
tent approximate. The obtained satisfactory agreement between the 
theoretical and the test results is a consequence of the smallness of 
the divergence angles of the diffusers investigated. Undoubtedly, with 
further increase in the divergence the effect of the secondary flows 
will show up to an increasingly greater extent and finally will lead to 
separation of the flow. The investigation of a separated flow in the 
diffuser is of great theoretical and practical interest. There does 
not as yet appear any way leading to the solution of this problem in 
the full sense of the word: that is, a solution which is self-contained 
with respect to the boundary conditions. In the solution of such a 
problem. further development of this method of considering the flow as 
emanating from a fictitious source may be of value. 

Noted here are a few of the most important conclusions drawn from 
the present investigation: 

1. The assumption made at the beginning of the investigation of 
the radial! ty of the flow holds true to a satisfactory degree for both 
diffusers investigated.. 

2. The assumption of the absolute character of the curve of mixing 
length for the straight pipe and diffusers and the representation of 
this curve by a cubical parabola give good agreement of the obtained 
velocity profiles with experiment. 

3. The resistance formula for the diffusers is identical with that 
for the pipe . 

4. The values of the absolute turbulence constants X and A in 
the resistance formula determined on the basis of tests on the straight 
pipe are applicable also to diffusers, and this confirms the absolute 
character of the constants. 

5. The approximate representation of the acutally nonradial flow 
in a diffuser by a radial flow originating from a fictitious source 
gives the necessary correction in computing the velocity drop along the 
axes of the 1° and 2° diffusers and in computing the pressure drop 
along the 1° diffuser. In computing the pressure drop for the 2° 
diffuser, less favorable results were obtained because sufficient ac- 
count was not taken of the magnitude connecting the normal component 
of the turbulent stresses with the gradient along the diffuser. This 
magnitude cannot as yet be theoretically obtained. 
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In conclusion, the writer wishes to express his deep appreciation 
to P. E. Euryatnikov for assisting in the tests and computations in 
connection with the present paper. 


Translation by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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Figure 2,- Friction distribution over a 

diverging, straight and converging pipe. 



Figure 3.- Correction coofficient £, taking account of the effect of 
the wall on the Pitot tube reading, plotted as a function 
of the distance from the wall. 



Variation of r^r 3 along axis of 1° diffuser for 
maximum discharge rate. 


Figure 8.- 
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Figure 9.~ Variation of t^r 3 along axis of 2° diffuser 
for maximum discharge rate. 
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Figure 11.- Variation of u m r^ along axis of 2° diffuser 
for various discharge rates. 
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Figure 12.- Pressure distribution along 1° diffuser, 
p = f(l/r 4 ). 
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Figure 14.- Comparison of theoretical with experimental velocity- 

profiles for 1° diffuser. 
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Figure 15.- Comparison of theoretical with experimental velocity- 

profiles for 2° diffuser. 
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Figure 16.- Comparison of the approxi- 
mation for the inertia 


integral I 



Id! 


with the results of graphical integra- 
tion from the theoretical velocity dis- 
tribution for the 1° diffuser. 


Figure 17.- Comparison of the approxi- 
mation for the inertia 



with the results of graphical infc ['ra- 
tion from the theoretical vein cl by 
distribution for the 2° diffuser. 
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1° diffuser. 
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Figure 



x for Ar m = - 2.27 m for 2° diffuser. 



Figure 29.— Dependence of S = — ^ on t m = 

— 2 P d C( r— Ar m 2l] 4 ] 

u m (r-Ar m ) fori 0 diffuser for Ar E = - 2.75 a. 



Figure 30.— Dependence of G = *- — on t = 

g q p d [( r- Ar m ) 4 ] “ 

n n ( r — ir E ) for 2° diffuser for Ar m = — 2.27 m. 










